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ABSTRACT 

A solution is presented for the spectrum of high-energy GRB photons confined to a quasi- 
thermal baryonic photosphere. The solution is valid in the steady-state limit assuming the 
region under consideration is optically thick to the continuously injected photons. It is shown 
that for a high luminosity photosphere, the non-thermal electrons resulting from 7-ray Comp- 
ton cooling lose their energy by upscattering the soft thermalised radiation. The resulting spec- 
tral modifications offer the possibility of diagnosing not only the burst comoving luminosity 
but also the baryon load of the ejecta. This model leads to a simple physical interpretation of 
X-ray rich bursts and anomalous low-energy slopes. 
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1 INTRODUCTION 

A quasi-thermal photosphere is expected in all GRB ejecta (Cavallo & Rees 1978; Goodman 1986; Shemi & Piran 1990), either due to pairs 
or to baryonic electrons, in addition to a possible non-thermal component from dissipation (shocks or reconnection). Preponderance of one 
or the other depends on whether the photosphere occurs inside or outside the saturation radius where the bulk Lorentz factor saturates to 
the dimensionless entropy of the outflow. Below the saturation radius, the photospheric luminosity dominates whereas above it the greater 
part of the energy is in kinetic form. In the latter case, the energy available in the from of radiation or pairs depends on the variable portion 
of the outflow and the efficiency of dissipation in shocks (Daigne & Mochkovitch 2002; Ramirez-Ruiz et al. 2002), or on the reconnection 
efficiency (Meszaros & Rees 2000). 

In this Letter we investigate the relationship between the quasi-thermal baryon-related photosphere in relativistic outflows, and the 
internal shock arising around this limiting radius. This photosphere is a source of soft thermal radiation, which may be observationally 
detectable in some GRB spectra (Murakami et al. 1991). In these models, dissipation happens whenever internal shocks develop in the 
ejecta (Rees & Meszaros 1994), which reconverts some fraction of the kinetic energy into radiation. If this dissipation takes place below the 
photosphere, the non-thermal shock luminosity could compensate for the adiabatic decrease of the post-saturation photospheric luminosity, 
so that the radiative efficiency in the outflow would be high. The relative roles of photosphere and shocks are reexamined in Section 2. A 
solution is presented in Section 3 for the spectrum of high energy 7-rays confined to a quasi-thermal baryonic photosphere and the possible 
non thermal spectral distortions in it. The relation of theses ideas to GRB phenomenology is outlined in Section 4. 



2 PHOTOSPHERES AND SHOCKS 

Consider a relativistic wind outflow expanding from some initial minimum radius ro = cSt — 10 7 ro,7 cm, where the wind baryon load 
M is parametrized by a dimensionless entropy 77 = Lq/Mc 2 . We then assume that the actual value of r/ (or Lq) is unsteady. The Lorentz 
factor saturates to T ~ r\ at a saturation radius r n /ro ~ r\ where the wind energy density, in radiation or in magnetic energy, drops below 
the baryon rest mass density in the comoving frame (Shemi & Piran 1990). The location of the baryonic photosphere, where tt = 1, due to 
electrons associated with baryons, is 

Mar _ ln i3 r -3 , n 
10 io,52??2 cm . (!) 



Anm v cT 2 



where 772 = 77/ 10 2 and L0.52 = Lo / 10 52 erg. The above equation holds provided that 77 is low enough that the wind has already reached its 
terminal Lorentz factor at r T . This requires (Meszaros & Rees 2000) 

/ Lqctt \ ^ ,„3 r i/4 -1/4 ,,,, 
77 ^={4^F ) - 10L oW- (2> 
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Figure 1. Schematic spacetime diagram in source frame coordinates of a relativistic outflow. The axes (logarithmic) are r versus t — (r/c), where t is time 
measured by a distant observer, and is zero when the burst is observed to start. In this plot, light rays are horizontal lines. The primary gamma-ray emission is 
assumed to continue, with a quasi-steady luminosity L w , for a time t w . If T fluctuates by a factor of ~ 2 around its mean value, relative motions within the 
outflowing material give rise to internal shocks. Decreasing 77 values lead to world lines further to the left. In case (a), r) < r) T and the dissipation occurs when 
the wind is optically thick. In case (b), with rj^ > r\ > r; T , ejecta collide in an optically thin region before reaching the contact discontinuity. The contact 
discontinuity and the forward shock are being decelerated because of the increasing amount of external matter being swept up (giving rise to a long-term 
afterglow; case [d]), so that they lag behind the light cone by an increasing amount Ar , whose increase with r is steeper than linear. This deceleration allows 
ejecta to catch up and pass through a reverse shock just inside the contact discontinuity (case [c]). 



If the value of 77 at the base increases by a factor > 2 over a timescale St, then the later ejecta will catch up and dissipate a significant 
fraction of their energy at some radius r L > r v given by 

r L ~ cStrf ~3x 10 14 <5to?72 cm. (3) 

The shocks cannot occur when r < r v since in this region both shells accelerate at the same rate T oc r and do not catch up. In order for a 
shock to occur above a photosphere which is in the coasting region 

r, > Vt « 50L$,ft^ 1/6 , (4) 

if one takes ro ~ cSt, where Sta = St/1 sec. If, on the other hand, St ~ 10~ 3 sec, then r] T ~ 200. 

The initial wind starts to decelerate when it has swept up ~ 77" 1 of its initial mass. For sufficiently high 77, the deceleration radius can 
formally become smaller than the collisional radius of equation This requires 

o 1 r,2 r 1/8 ,1/8 -1/8 r ,-3/8 , c , 

77 > r?d « 8 x 10 L ' i52 tJ A n St . (5) 

This deceleration allows slower ejecta to catch up, replenishing and re-energising the reverse shock and boosting the momentum in the blast 
wave. 

Fig. 1 shows the schematic world-lines of a relativistic outflow with a range of Lorentz factors. We identify three types of contributions 
to the observed time history, each with a different character. For a relatively low Lorentz factor 77 < r) T , as in curve (a) of Fig. 1, the radius 
r l < rv, and the dissipation occurs before the wind is optically thin. Below the baryonic photosphere, shocks would occur at high optical 
depths. Section 3 investigates the relationship between the quasi-thermal baryon-related photosphere in relativistic outflows and the internal 
shocks arising near this limiting region. For a larger Lorentz factor 77 > r/d, corresponding to curve (c) of Fig. 1, the ejecta would expand 
freely until the contact discontinuity had been decelerated by sweeping up external material. It would then crash into the reverse shock, 
thermalising its energy and boosting the power of the afterglow. The impact of such collisions on the prompt emission have been studied in 
Ramirez-Ruiz, Merloni & Rees (2001). In curve (b), with intermediate 77, deceleration occurs at radii r& > r L , and dissipation takes place 
when the wind is optically thin (i.e. when it is most effective). 



2.1 Photospheric and Internal Shock Luminosity 

The lab-frame baryonic photospheric luminosity L T and dimensionless temperature T evolve as (L T /Lo) = (G r /6o) = (?V /r v )~ 2 ^ 3 , 
where Oo = kTo/m e c 2 ~ 2L\^ 2 r is the initial temperature at ro. Using the expression for r T , this can be written as (L T /Lo) = 
(rj/v*) 8 ^ 3 for?7 < 77* and 1 otherwise. 

The internal shocks in the wind can dissipate a fraction of the terminal kinetic energy luminosity Lrj above the saturation radius r v , 
L L = e e e L Lo ~ 10 _1 (3e e )(3e t )io , where e e e t is a bolometric radiative efficiency when the cooling timescale is shorter than the dynamical 
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time. If the Poynting flux provides a fraction a of the total luminosity Lq at the base of the wind (at ro), the magnetic field there is 
B ~ 10 10 a 1/2 L^ 5 2 1 <5*Q 1 G. The comoving magnetic field at r L is B L = Bo(r /r v ) 2 (r v /r L ) ~ 10 4 a 1/2 L^^St' 1 ^ 3 G. If the electrons 
are accelerated in the dissipation shocks to a Lorentz factor 7 = 10 3 73, the ratio of the synchrotron cooling time to the dynamic expansion 
time in the comoving frame is (tsy/tadi)t ~ 5 x 10~ 3 a _1 L ~g 1 7^ 1 (5fo?72 so a ver Y high radiative efficiency is ensured even for St as high 
as seconds. It is therefore clear that a magnetic field can ensure efficient cooling even if it is not strong enough to be dynamically significant 
(i.e. even for a <C 1). 

For a given Lo, an individual burst is characterized by an average 77. For 77 > 77*, the photosphere arises in the accelerating region, and 
in this region shocks are not possible. For r\ T < 77 < 77,, we have L T > L L , so that the baryonic photospheric component dominates the 
non-thermal internal shock component in a bolometric sense. Low values of 77 lead to further out, weaker photospheres and, at the same time, 
relative stronger shocks occurring near the photosphere. For L T < L L , the thermal photospheric peak will tend to blend with the synchrotron 
peak, resembling the canonical non thermal GRB spectrum. For lower 77 < rj T , the thermal peak is even softer with L r <Lj, while shocks 
occur closer in and produce harder synchrotron peaks (Ramirez-Ruiz & Lloyd-Ronning 2002). The corresponding shocks occur below (or 
near) the baryonic photosphere where the scattering depth of the shells is larger (Panaitescu et al. 1999; Ramirez-Ruiz & Lloyd-Ronning 
2002, Kobayashi et al. 2002). Shocks which occur inside the photosphere may also induce Alfven waves and these waves can be efficiently 
damped for r < r T (i.e. tt > 1). 



3 THE EFFECTS OF A REPROCESSING PHOTOSPHERE ON THE GRB SPECTRA 

We consider here the problem of Compton downscattering of X-rays and 7-rays confined to an optically thick baryonic photosphere. The 
electrons are assumed to have zero temperature, which is appropriate as long as the photon energies are much larger than T . For simplicity 
we consider here only shocks producing a comoving spectrum with hv' < 511 keV so that photon-photon pair production is negligible. 
In this case, the pair photosphere is not sufficient to alter significantly the spectrum of the baryonic photosphere. In the following we shall 
assume the acceleration of electrons in the shocks to be impulsive, and therefore it has to take place in a very limited volume of the interacting 
shell (i.e the emitting particles do not provide an additional heating source). 

We assume that GRB photons are continuously and uniformly injected at a rate 0i(e), where e is the dimensionless photon energy in 
units of m e c 2 . For future use we take a power-law photon injection of the form 4>i(e) = ^/e~' 3 , for e < et. Under the assumption that 
photosphere is optically thick, a steady photon density distribution <f>(e) can be found after specifying the photon production rate and the 
electron density n e of the region under consideration. This distribution is the solution to a kinetic equation where the main energy-loss 
mechanism is Compton scattering. <j>i(e) satisfies an equation of the form (Arons 1971) 

= &(e) + f " efeV(e>(e', e) - cj>{e) [ de'a(e',e), (6) 

where the three terms correspond to the injection of primary photons, Compton downscattering to e from higher energies, and Compton 
scattering out of the energy e to lower energies, respectively. <r(e', e) in equation J6j represents the probability per unit time that a photon 
with energy e will Compton scatter from e' to e (Jauch & Rohrlich 1980) 



e' e /l 

7 + 7- 



+ 



(7) 



The integral equation |6| with the scattering rate can be solved numerically by iteration for any e, where the iterative steps correspond 
to subsequent orders of Compton scattering. In Fig.|2|are plotted representative spectra for a power-law injection of high-energy photons in 
the conventional coordinates e and e 2 <^(e). The dotted curves give the density of the unscattered photons, (f> u , produced at the power law rate. 
The solid curves, on the other hand, give the total photon density which includes both unscattered and scattered photons. Here cf> u is given by 
<^u(e) = ^i(e)/[tief(e)c]. The calculations shown in Fig. [^approximate the electrons in the medium to have zero temperature and assume 
that Compton scattering is the main interaction process for 7-rays. The results shown in Fig.[3]expand this calculation to include the role of 
secondary electrons (produced by Compton downscattering) under the assumption that they lose their energy mainly by upscattering the soft 
fhermalised radiation (see Section 3.1). 

The right panel of Fig.|2|shows the evolution of the photon density distribution for subsequent orders of Compton scattering. The black 
solid curve gives the total photon density distribution which includes both unscattered and scattered photons. The grey solid line corresponds 
to the density of unscattered photons, and the grey dashed curve to the sum of unscattered and singly scattered photons. The density of 
photons scattered once and more is therefore given by the difference between the black solid and grey dashed curves. 

When e <C 1, it is straightforward to present a simplified discussion of the problem. Analytical solutions are derived here largely to 
illustrate what may not be obvious from the numerical treatment. The last integral in equation ^6) can be rewritten as 

de'a(e, e) = n e ca(e), (8) 

/(l+2e) 

where cr(e) is the Klein-Nishina cross section. In the case of e <C 1, <r(e) ~ <jt- The Fokker-Plank equation can then be derived by expanding 
^6) in powers of the small quantity e — e' (Ross et al. 1978) 
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Figure 2. The photon density distribution from a power-law injection of 7-rays. The normalizations of <f>(e) correspond to \& /(n e a'T'c)=l. Left Panel: Power- 
law injection of primary photons with index (3 and £{,=5 (solid lines). The dotted curves correspond to the density of unscattered photons. Right Panel: The 
evolution of the photon density distribution in the presence of a power-law injection of primary photons with index /3 = 1.5 and e(,=10. The black solid line 
gives the total photon density distribution. The grey solid line corresponds to the density of unscattered photons, and the grey dashed curve to the sum of 
unscattered and singly scattered photons. 
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where the first term describes the photon energy dispersion, whereas the second characterizes the systematic changes in e. Equation {5) differs 
from the equation of Kompaneets (1957) in that there is an additional dispersion term. For a power-law injection, the solution to equation |9j 
reads 



n e OTc{[5 ■ 



,1-/3 



7e{/3 + 3} 
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Here we have kept only the first two powers of e. Provided that e\, > 1 and (3 > 1, from equation i 101 it follows that 
illustrated in Fig.|2| most of the luminosity at low energies results from 7-ray downscattering. 



(10) 



As 



3.1 The role and fate of the secondary, energetic electrons 

So far we have approximated the electrons in the medium to have zero temperature and have assumed that Compton scattering is the main 
interaction process for 7-rays. Compton cooling of 7-rays, however, produces a population of secondary energetic electrons. The fate of 
these electrons depends on the composition, ionization state and magnetic field present in the medium, as well as on the energy density of the 
background radiation. When the energy density in the (black body) photospheric photons is much larger than the energy in scattered photons 
(i.e. L T > L L ), the secondary electrons will lose energy mostly in collisions with the photospheric photons and form a secondary photon 
spectrum. This is the situation considered in this section. 

To first approximation, the scattering medium can be considered as a cold plasma if the time between successive scatterings, t T ~ 
em e c/ '{UlO't), is much longer than the characteristic cooling time of the secondary electrons, t 1 ~ jm e c/(-yU T crT)- Here U L and U T are 
the energy densities of the continuously injected photons and the photospheric radiation, respectively. The secondary electrons created by 
downscattering Compton cooled before the next generation of hot electrons is created if U T > jU L /(e / y). Therefore, the bulk of the injected 
radiation is likely to be scattered by cooled particles when L T > L L . In what follows we assume the scattering medium to be "cold" since it 
greatly reduces the calculations. A detailed discussion of the relevance of a hot scattering medium, its steady state solutions for various input 
spectra, and its range of applicability (i.e L T < L L ) will be presented in future work. 

The injected photons repeatedly downscatter to form a steady state distribution. During this process, a population of secondary, energetic 
electrons is injected. Down-scattering of 7-rays on a cold medium produces secondary electrons at a rate 

n 7 ( 7 ) = jf * de0(e)a(e,e-j + l), (11) 

where 7 is the electron Lorentz factor and <j> is the steady state photon density distribution. Here e m = §[7 — 1 + {(7 — l) 2 + 2(7 — l)} 1 ^ 2 ] 
is the minimum energy required for a photon to up-scatter an electron from rest to the energy 7. 

The maximum energy of the secondary electrons is then given by 75 = 1 + 2ef /(l + 2eb). These secondary electrons could lose their 
energy mainly through synchrotron or Compton energy-loss mechanisms. One or both of these mechanisms may be present, depending on 
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Figure 3. Illustration of the effect of Compton scattering of photospheric blackbody photons at 7V =0.5 keV by secondary electrons. The role of secondary 
electrons is included under the assumption that they lose their energy mainly by upscattering the soft thermalised radiation and that they Compton cooled before 
the next generation of hot electrons is created. The production rate of secondary photons resulting from scatterings of the secondary energetic electrons with 
the black body photons is calculated and subsequently added to the total rate of photon injection. This procedure is repeated until convergemce is achieved. 
The dark solid line gives the spectra without the effects of blackbody up-scattering (Fig.|5J. The grey solid line, on the other hand, illustrates the spectral 
distortions arising from the inclusion of photospheric photons while the dotted line gives the photon distribution after only the first iteration procedure. The 
primary photons are injected with a power-law distribution with fj = 1.5 and ej, = 10. The normalization of <j>(e) correspond to (n e a-yc)=l. 
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the bulk Lorentz factor, the isotropic equivalent total energy of the burst, and the efficiency of dissipation in shocks. In the previous section 
we discussed the criteria for the non-thermal component to dominate over, or be dominated by, the photospheric thermal component. A 
photosphere should be prominent in the lowest baryon load cases. The black body photons can then act as seeds for scattering to higher 
energies, if there is substantial amount of energy in scattering centers. Alfven waves generated by magnetic field reconnection or MHD 
turbulence can act as such centers. Repeated scattering on the Alfven waves acts in the same way as Comptonisation off hot electrons. 

Electrons resulting from down-scattering will therefore lose energy mostly in collisions with the photospheric emission when the energy 
density in the photospheric blackbody exceeds the scattered one. The steady state electron distribution (Blumenthal & Gould 1970) is given 
by 

717(7) = 4 7 roo , (12) 

37 2 (T T c J Q de e ee0e(ee) 

where 0e(ee) is the blackbody photon density distribution. 

The production rate of secondary high-energy photons (i.e. 6 > O r ) produced by the upscatter of blackbody photons by the secondary 
electrons is given by 

rib 

4> s {e)= I ^7^(7)0-0(7,7 - e)> (13) 

Jl + e 

where <j@("/, 7') is the rate of scatterings of electrons from 7 to 7' on blackbody photons. Zdziarski (1988) has given an explicit expression 
for ae (7, 7') • The assumption of negligible photon-photon pair production guarantees that collisions between secondary electrons and black 
body photons take place in the e<l limit. The total rate of photon injection is then given by (f> s (e) + <pi (e) . This photon rate results in a new 
steady state photon distribution J6j, and, in turn, a new secondary electron production rate 1111 . The production rate of secondary photons 
resulting from scatterings of the new secondary energetic electrons with the black body photons is calculated and subsequently added to the 
total rate of photon injection. This procedure is repeated until convergence is achieved (i.e. (j)™ + <j)" for n — 1,2, 3, ...). For each step the 
integral equation {6j is solved with the scattering probability {7) numerically by iteration (see e.g. right panel of Fig.|2j- 

Fig-EJillustrates the effect of Compton up-scattering of photospheric blackbody photons at G T ~ 10 _3 m e c 2 by secondary electrons. 
The dotted curve shows the steady state photon distribution with the inclusion of the soft black-body photons while the black dashed line 
gives the photon distribution after only the first iteration procedure. The photons from photospheric upscattering dominate the spectrum 
at e <C 1. The highest energy (upscattered) photospheric photons are produced by relativistic electrons from Compton downscattering of 
the primary photons. At lower energies, on the other hand, the upscattered photospheric photons are due to electrons arising from multiple 
scattered photons. The importance of the upscattered photospheric photons at e < 0.1 depends mainly on the low-energy tail distribution of 
the primary photons and thus its effect on the photon spectra increases with subsequent iterations. This emission could be responsible for 
the X-ray excess observed in some GRBs. In contrast, at higher energies, the injected photon distribution (i.e. <fi™ + 0") is not significantly 
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alter by subsequent iterations and the role of upscattered photospheric photons to the high-energy spectra can be estimated with only a few 
iterations. 

One does not expect this scattered photospheric component to be a pure Wien spectrum, since in order to change the shape of a soft 
spectrum and make it into a dilute black body spectrum would require that each photon be scattered m e c /hv times by a population of 
electrons in thermal equilibrium. This requires an optical depth much larger than unity, which is generally not the case anywhere near the 
photosphere. The scattering depth per shock due to pairs is also unlikely to be much larger, because the scattering and the pair-formation 
cross sections are comparable, and unless dissipation and pair formation occurs uniformly throughout the entire volume, down-scattering of 
photons above the pair threshold rapidly leads to self-shielding (Meszaros et al. 2002). 



4 DISCUSSION 

The standard internal shock model of GRB is generally assumed to produce its observed non thermal radiation by synchrotron (or possibly 
inverse Compton) process. Here, in addition to synchrotron we have also considered in more detail the role of the outflow photosphere and of 
possible non thermal distortions in it. A strong baryonic photosphere component should be present at the beginning of bursts. However, the 
farther beyond the coasting radius the photosphere occurs, the weaker its energy fraction is relative to the primary injected photons arising 
from internal shocks, because its energy drops as r~ 3 / 2 . 

The problem of reprocessing X-rays and 7-rays by Compton scattering could be the key to understanding the X-ray excess above the 
power-law extrapolation from higher energies observed in a non-negligible fraction of bursts (Kippen et al. 2001; Heise et al. 2001; Amati 
et al. 2002; Lloyd-Ronning & Ramirez-Ruiz 2002; Sakamoto et al. 2005), although smoothing and softening could be stronger where there 
is substantial pair formation (Pilla & Loeb 1998; Ghisellini & Celotti 1999; Meszaros et al. 2002; Pe'er & Waxman 2004). If there were a 
power-law non-thermal spectrum in which significant fraction of the resulting photons are above 1 MeV in the comoving frame, then pair 
production will change the situation. Photons with a few MeV energies in the comoving frame will then be converted into pairs with very 
modest 7 provided that the compactness parameter is more than unity. The pairs will generally have relativistic energies, and will themselves 
participate in the synchrotron and inverse Compton emission (Pe'er & Waxman 2004). In reality the actual time dependence for an unsteady 
outflow leading to shocks, pair formation, and Comptonisation could be more complicated. 

If both a photospheric and a shock component are detected, one would expect the thermal photospheric luminosity (and its non thermal 
part, if present) to vary on similar timescales as the non thermal synchrotron component (unless the shock efficiency is radius dependent, 
or unless one or both are beyond r — rorf , in which case r/(cq 2 ) imposes a lower limit on the corresponding variability timescale). The 
luminosity in a given band (e.g. Swift) probably varies differently, since the thermal peak energy is tx L 1 / 4 and falls off steeply, while the 
synchrotron peak energy varies L A/2 and falls off more slowly. A preferred low-energy break (Fig. 3) could be attributed to a scattered 
photospheric component, provided the baryon loads are low or the outflow variability timescales are large (Fig. 1). This still requires a 
relatively strong shock synchrotron component, or possibly Alfven wave Comptonisation in the photosphere (e.g. Thompson 1994), to 
explain the high-energy power-law spectra. It would also imply a pronounced upward change of slope above the X-ray excess, from the 
thermal peak to a flatter power law in all bursts where a low break is observed. Further data on X-ray and 7-ray spectral features during the 
burst (as opposed to the afterglow), will surely offer important clues to the nature of the bulk flow and the macroscopic source of energy 
driving the microscopic processes of particle acceleration and cooling. 
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